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Abstract 

In this article we present the cut Fock space approach to the D = 
d + 1 = 2, Supersymmetric Yang-Mills Quantum Mechanics (SYMQM). 
We start by briefly introducing the main features of the framework. We 
concentrate on those properties of the method which make it a convenient 
set up not only for numerical calculations but also for analytic computa- 
tions. In the main part of the article a sample of results are discussed, 
namely, analytic and numerical analysis of the D — 2, SYMQM systems 
with SU(2) and SU(3) gauge symmetry. 

1 Introduction 

Supersymmetric Yang-Mills Quantum Mechanics (SYMQM) turned out to be 
not only a class of supersymmetric systems possessing some interesting phys- 
ical features, but also to play an important role in many areas of theoretical 
physics. Among the problems where SYMQM are relevant, the two most no- 
table examples are: their relation to a particular limit of M-theory pQ, and 
the description of regularized dynamics of relativistic quantum membranes and 
supermembranes [5] [5] . Hence, an efficient way of investigating the spectra of 
SYMQM with various gauge groups and defined in spaces with different dimen- 
sionality would be of importance. 

In short, D = 2, SYMQM are supersymmetric, D — 2 dimensional, Yang- 
Mills quantum field theories reduced to one point in space. Therefore, the 
original local gauge symmetry is transformed into a global symmetry of the 
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quantum mechanical system. The difficulty of solving these systems even in the 
simplest cases comes from the singlet constraint which is the remnant of the 
Gauss law. 

The cut Fock space approach [3] was proposed few years ago as a nonper- 
turbative way of investigating numerically the SYMQM systems. The energy 
eigenstates are constructed as linear combinations of physical basis states and 
thus are gauge invariant by construction. The approach allowed the calculation 
of low lying eigenenergies and their eigenstates through a numerical diagonal- 
ization of the Hamiltonian matrix. Its main application was the study of D = 2, 
SYMQM [5] [5] [7] , but it was also generalized to the D — 4 case [5] ■ In principle 
the method can be used in numerical investigations of systems with any gauge 
group and in any dimension. Nevertheless, we would like to stress in this article 
that it is also convenient for analytic treatment. The cut Fock space approach 
provides a systematic way of controlling the Fock basis in terms of which the 
energy eigenstates can be decomposed. As an example, we will present such 
decomposition of the solutions for two simple systems. 

The paper is composed as follows. We start by briefly describing the cut 
Fock space approach, concentrating on the construction of the Fock basis. Then, 
we present the numerical algorithm together with the numerical spectra of the 
D = 2, SYMQM with SU(3) gauge group in the sectors with n F = and 
np = 2 fcrmionic quanta. In the third part of the paper, we derive analytically 
the bosonic solutions of the D = 2, SYMQM Hamiltonians with SU(2) and 
SU(3) gauge groups. We conclude by indicating possible directions of further 
studies. 

2 The framework 

For the reasons of simplicity the framework will be presented in the context of 
D = 2, supersymmctric Yang-Mills quantum mechanics with the SU(N) gauge 
group. Nevertheless, the method is more flexible and systems with other gauge 
groups as well as in higher dimensional spaces can be studied. Particularly, 
creation and annihilation operators introduced in the following subsection can 
be labeled by additional spatial indices and transform in some representation of 
the SO(d) group. 

2.1 Basic degrees of freedom 

A two-dimensional SYMQM system is described [9] by a bosonic variable 4>a and 
a complex fermion A^, where A is a color index. Being remnants of the gauge 
field in the original field theory, the bosonic and fermionic variables transform 
in the adjoint representation of the SU(N) group. Thus, the system contains 
N 2 — 1 bosonic and N 2 — 1 fermionic degrees of freedom. 

The characteristic feature of SYMQM is that its Hilbert space is composed 
of states invariant under the SU(N) group. We would like to incorporate this 
constraint in the approach from the beginning. Hence, we introduce a matrix 
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notation, in which any singlet can be written in terms of traces of appropriate 
matrices HOl. We define 



AT-1 N'-l 

k] = 



A=l A=l 

N 2 -l N 2 -l 



A=l A=l 

where are the generators of the SU(N) group in the fundamental represen- 
tation, i,j — 1,...,N. Now, all operators become operator valued matrices. 
Particularly, the gauge-invariant occupation number operators can be defined 

as 

N 2 -l N 2 -l 

tr (ata) = Yl a A"A, tr (p f) = £ f\f A , (1) 

A=l A=l 

where we have introduced standard commuting creation and annihilation op- 
erators a A ,a A defined by 4> A = -^(a A + a A ) and ir A — ~^( a A ~ a A)- The 
sum over the adjoint indices ensures the invariance of these operators under the 
SU(N) transformations. 

In the following we will consider the D = 2, SYMQM systems with SU(2) 
and SU(3) gauge groups. Their Hamiltonians read respectively [5], 

1 3 1 

H = -n Air a = tr {a) a) + -- - (tr (aV) - tr (aa)), (2) 

H = ^ir A n A = tr (a f a) + 2 - i(tr (aV) - tr (aa)). (3) 
2.2 Fock basis 

The fundamental part of the approach is a systematic and recursive construc- 
tion of the Fock basis. Fock states are eigenstates of some occupation number 
operators, and in the case of SYMQM models, we choose them to be eigenstates 
of the gauge- invariant occupation number operators eq.(P). 

It is convenient when the Fock basis states can be labeled by as many quan- 
tum numbers conserved by the Hamiltonian as possible. Hence, since for most 
of systems the fermionic occupation number is conserved, one usually constructs 
the Fock basis independently in each subspace of the Hilbert space with a defi- 
nite fermionic occupation number. As far as the bosonic occupation number is 
concerned, it is in general not conserved. Nevertheless, we will further divide 
the fermionic sectors into subspaces with a given number of bosonic quanta, in 
order to facilitate the recursive approach. 

In the following we introduce the concept of bosonic elementary bricks, which 
are necessary to obtain the Fock basis in the bosonic sector. Then, we will 
proceed in full analogy with the fermionic sectors. 



3 



2.2.1 Bosonic elementary bricks 



We define the set of bosonic elementary bricks as the set of iV linearly indepen- 
dent single traces of bosonic creation operators. Traces with more than N — 1 
operators can be reduced by the Cayley-Hamilton theorem. Table Q] contains 
examples of such sets for N = 2, N = 3 and N = 4. 



517(2) 


517(3) 


517(4) 


tr (aV) 


tr (at at) 
tr (aWa^) 


tr (atat) 
tr (a t a t a t ) 
tr (atatatat) 



Table 1: Elementary bosonic bricks for SU(2),SU(3) and 5*7(4). 



Let us consider the set of stated 

{(at 2 )* 2 (at 3 )* 3 ...(a tAr ) feN |0)} 



5Zj = 2 3&3 ~ 



= \{n B }) 



(4) 



composed of the products of powers of elementary bosonic bricks acting on the 
Fock vacuum. One can show [7] that it spans the subspace of the Hilbert space 
with tlb bosonic quanta. In eq. ([!]), we introduced a generalized notation in 
which |{ns}) is a vector of all states with ub quanta. 

Suppose that we have constructed such basis up to sectors containing less 
than ub bosonic quanta. The Fock basis in the sector with ub bosonic quanta 
can be build as the sum of all states obtained by the action of appropriate bricks 
on the already generated Fock basis states. Using the generalized notation this 
can be written in a compact form as 



|{n B })=]T(at fc )|{n B -fc}). 



k=2 



Note that the same state may appear in several copies, differing in the order of 
successive bricks used to build it. The basis is obtained once this redundancy 
is removed and the remaining states orthonormalized. 



2.2.2 Fermionic bricks 

In order to obtain the Fock basis in the fermionic sectors we must define bricks 
which contain fermionic creation operators. The set of elementary fermionic 
bricks can be defined in full analogy to the set of elementary bosonic bricks. We, 
thus, consider all single traces with np fermionic creation operators, which can- 
not be further reduced by the Cayley-Hamilton theorem. Subsequently, such set 
of elementary fermionic bricks must be enlarged by operators, which are prod- 
ucts of fermionic elementary bricks with smaller number of fermionic quanta and 

1 We adopt here the notation in which (X) designs tr (X). We will use this notation only 
when it is self-evident. 
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F = 1 


_F = 2 


F = 3 


F = 4 


(/ T <* T ) 
(/t t t) 


(/t/tot t) 
(ft a t a t/t a t) 
(/t a t)(/t a t a t) 


(/W) 

(/WW) 
(/Wx/W) 

(/tflt/t/t a t tt t) 

(/t a t)(/t/t a tat) 
(/t a t a t)(/t/t a t) 
(/t a t a t)(/t / t a t a t) 


(/WW) 
(/ f at)(/t/t/t) 

WWW) 

(/t a t a t )(/ t / t / t) 
(/+ot)( t/t/t/t) 

(/wx/vw) 

(/t a tat)(/t/t/t a t) 

(/ t / t at)(/t/t a t a t) 

(/W)(/W)(/W) 
(/t/tot)(/t t/t a t t) 



Table 2: Fermionic SU{3) bricks. 

contain np fermionic creation operators in total. The inclusion of these com- 
posite operators ensures that all possible invariant contractions of np fermionic 
and some bosonic creation operators are taken into account. The enlarged set 
is called the set of composite fermionic bricks. Its elements will be denoted by 
C' (n,B,riF,a), where the first argument describes the number of bosonic cre- 
ation operators contained in C', and the last argument is an additional label 
needed in the cases where tib and uf are not enough to distinguish different 
operators. Table [2] presents the set of composite fermionic bricks for the SU(3) 
gauge group. 

In analogy to the bosonic case, we can define the set of states, 
{cHn,n F ,a)(a^(a^...{a^) k -\0)} =\{n B ,n F }) (5) 

which after orthonormalization will give the basis in the subspace of Hilbert 
space with hb and np bosonic and fermionic quanta, respectively. 

2.2.3 Correctness of the Fock basis 

One can show that the sets of states eqs. (|4]) and (O provide indeed a good 
basis of the Hilbert space. On one hand, its completeness follows from the fact 
that the states constructed with powers of bosonic and fermionic bricks repre- 
sent the most general contractions of invariant tensors with creation operators. 
On the other hand, after removing the trivially linear dependent states, the 
linear independence of the remaining ones can be check explicitly by calculating 
the Gram's matrix. Obviously, the rank of the Gram's matrix corresponds to 
the number of linearly independent states in a sector with given ng and np. 
Fortunately, there exists also an independent way of calculating this number 
PTj . It exploits the orthogonality of the characters of the SU(N) group, and 
can be used as a crosscheck that the Fock basis obtained through the recursive 
construction spans correctly the physical Hilbert space of SYMQM. 
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2.3 Extraction of approximate eigenenergies 
and eigenstates 

Once the Hamiltonian operator is expressed as an operator function of creation 
and annihilation operators, its action is straightforward in the Fock basis. How- 
ever, the numerical analysis requires an additional step, namely the introduction 
of a cut-off N cut (see figure Q]) on the countably infinite Fock basis. There are 
many ways to introduce such a cut-off depending on the symmetries of the sys- 
tem. A practical cut-off is a limit on the total number of quanta contained in the 
Fock basis states. Once the cut Hamiltonian matrix is obtained, its eigenvalues 
correspond to an approximation of the eigenenergies of the quantum system, 
and its eigenvectors to the eigenstates. Finally, calculations with several in- 
creasing N cut have to be performed and the physical results extracted from the 
limit of infinite cut-off. The properties of such a procedure were analyzed in 
[12] [13] [14] [15] , where a different behavior of the eigenenergies corresponding to 
localized states and those corresponding to nonlocalized states was observed. 
Therefore, the method offers a tool to distinguish these two types of states, and 
we will indeed exploit this possibility when discussing the numerical results in 
subsection 13.21 
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Figure 1: Physical interpretation of the introduced cut-off. At finite cut-off 
the position and momentum operators have discrete spectra. As an example 
eigenvalues of the pair of operators 4>i and </>2 as well as 7Ti and 1T2 are shown for 
different cut-offs (N cut = 5, 11 for the left and right figures). With increasing 
cut-off the points become more and more dense and tend to cover a bigger area 
on the plane. In the limit of infinite cut-off the operators have a continuum 
spectrum and the whole plane is covered. 
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3 Numerical results 



In this section we briefly describe a recursive algorithm which can used to ef- 
ficiently evaluate the matrix elements of the Hamiltonian operator. Then, we 
present the spectra of D = 2, SYMQM system with SU(3) gauge group in the 
sectors with np — and np = 2. 

3.1 Recursive numerical approach 

The recursive algorithm is based on relations connecting the desired matrix 
element of an operator to simpler matrix elements of some other operators, 
which have been already evaluated at the earlier stage of calculations. We 
will not describe here the full algorithm in details (for a full presentation of 
the approach see |16j). Instead, we will concentrate exclusively on the bosonic 
sector and give one example of such recursive relation. 

Let us now assume that we want to evaluate the matrix elements of the 
normally-ordered operator 0(n B ) between states containing ns and n' B quanta, 
( { n B } 1 ( n B ) I { n B } ) ■ The argument of O means that the difference between the 
number of creation operators and the number of annihilation operators is n B . 
Therefore, the matrix elements ({^bII^K^s)!! 7 ^}) wm " ^ e nonzero only when 
n' B = n B + ns ■ The strategy to evaluate a matrix element of O is to drag O 
over the operators constituting the state so that O annihilates the Fock 

vacuum. We thus have 

(K}|0(ng)l{»B»=(<H}l[0(ng),(o t, ')]l{»fl-p}> 

+ ({n'B}\(a^)0(n°)\{nB - p})) ■ R(n B ), (6) 

where the matrix R{ns) is obtained from the matrix of scalar products S(ub) = 
and is used to remove redundant basis vectors and orthonormal- 
ize the remaining ones. Thus, R(ub) satisfies 

R T (n B )S(n B )R(n B ) = I, 

where I is the unity matrix which dimension is equal to the size of the subspace 
of the Hilbert space with ns bosonic quanta. Note that we have expressed the 
desired matrix element through eq. ^ in terms of matrix elements of operators 
between states with lower number of quanta, which should have been already 
evaluated during some previous calculations. Hence, applying successively such 
relation one can evaluate {{n' B } \0(n%)\{n B }). O ne can similarly organize the 
calculations in the fcrmionic sectors. 

3.2 SU(3) model 

As an application of the above algorithm we present the spectra of the sys- 
tem given by the Hamiltonian eq. ([3]) . Figures [2] and [3] contain the results in 
the sectors with np = and np — 2, respectively. We notice that in both 
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Figure 2: Dependence of the eigenenergies on the cut-off for the SU(3) model 
in the rip = sector. 



cases all eigenenergies fall down to zero. Thus, we can conclude that the corre- 
sponding eigenstates are nonlocalized. It is a feature of the approach that the 
eigenvalues corresponding to such states do not converge with the increasing 
cut-off. Such behavior results from approximating a plane wave by a finite set 
of localized harmonic oscillator eigenstates and its dependence on N cut can be 
parameterized in a power like manner. A precise procedure for obtaining the 
correct infinite cut-off limit of the energy was described in [12] [13] [15] . In par- 
ticular, it was shown that the hyperbolic fall off contains information about the 
dispersion relation. Note that for a finite N cut there can be only a finite number 
of eigenenergies. With increasing N cut more and more new eigenvalues should 
appear on plots figs. [2] and [3] in the limit giving the continuum spectrum of 
a free system. In the next section we will present analytic calculations which 
enable one to completely reconstruct the above spectra, and thus provide a full 
understanding of their features. 



4 Analytic solutions 

Our procedure to derive exact solutions for the D = 2, SYMQM models consists 
in three steps. We decompose a general state \E) from the cut Hilbert space 
in the Fock basis. In this step a convenient parametrization of the Fock basis 
is necessary. Then, we translate the requirement that \E) is an eigenstate of 
the appropriate Hamiltonian into a recurrence relation on the decomposition 
coefficients. The last step is the solution of this recurrence relation at finite 
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cut-off and the investigation of the infinite cut-off limit of these solutions. 

Before we deal with the Hamiltonian ([3]), we will analyze a simpler model 
with the SU(2) symmetry group, ([2]). The solution of the SU (2) model was first 
found by Claudson and Halpern [9 . It can be obtained from the generalized 
solutions of Samuel [T7] and was also recently rederived by Trzetrzelewski [7] by 
algebraic methods. This model is particularly simple because the SU{2) group 
is isomorphic with the SO(3) group and can be parameterized by spherical 
coordinates. 

4.1 £77(2) model 
4.1.1 Recurrence relation 

At cut-off N cut the general state from the bosonic sector can be decomposed as 

N cut 

\E) = Y / a j (E) tr (a^O). 
j=o 

The eigenequation 

H\E) = E\E), 
yields the recurrence relation for the aj(E) coefficients, 

aj -i{E) - (2j + | - 2E) aj (E) + (j + + \)a j+l {E) = 0, (7) 
which can be solve analytically as will be discussed in the following subsections. 
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4.1.2 Finite cut-off solutions 

One can show [T5] that eq.® admits N cut + 1 solutions. The possible eigenen- 
ergies are given by the zeros of an appropriate associated Laguerre polynomial, 

4 cut+ i(2£)=0. 

For each E satisfying the above equation there exist an eigenstate with N cut + 1 
decomposition coefficients are given by 

a 3 {E) = a T(^)L](2E) 0<j< N cuU 
where ao is some arbitrary constant. 

4.1.3 Infinite cut-off solutions 

In the infinite cut-off limit, N cut — > oo, the set of possible eigenenergies is given 
by the whole real positive axis. This reflects the fact that the physical spectrum 
is continuous. Hence, for any real number E, there exists an eigenstate of H 
which decomposition coefficients are given by 

aj (E) = aoT(~)L](2E) j > 0. 
Therefore, the exact eigenstate can be written as, 

„ oo i 

\E)=a o n-)Y,L](2E) tr (a^O). 

4.1.4 Reconstruction of wavefunctions 

In the case of the SU(2) model one can explicitly check the correctness of the 
above solutions. Let us denote by ip n {r) the wavefuntion in the position repre- 
sentation of the n-th basis state, r is the radial variable which parameterizes the 
SU(2) group manifold. Since we consider only SU(2) invariant states, VVtM do 
not depend on angular variables. The n-th basis state being the eigenstates of 

2 Laguerre polynomials C„(x) are denned as the solutions of the differential equation 
xy" + (a + 1 — x)y' + ny = 0, 

and the orthogonality relation 

roo 

/ CZ l (x)C™(x)x a e- x dx = S mn . 
Jo 

L"(x) denotes the Sonine polynomials related to Laguerre polynomials by 

r(a + n + lj 
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the gauge-invariant particle number operator tr (aJa), its wavefunction satisfies 
the following equation 

1 / d 2 2d \ , , . , . , 

Such equation can be solved yielding 

Mr) = a(n)e~^(~l) n 2 2n+1 nld (r 2 ) + pe~$ -i-Fi(-^±l J, r 2 ), 

r 11 

where a(n) is some constant depending on n and ii 7 'i(a, 6, z 2 ) is the Rum- 
mer's function of the first kind. We are only interested in normalizable solu- 
tions, therefore since the function -iFi(a,b, z 2 ) is singular at z = 0, we set 
(3 = 0. Thus, V , n(r) turns out to be the wavefunction of a three dimensional 
harmonic oscillator carrying zero angular momentum. In order to find explicitly 
the Claudson-Halpern solutions, we write the eigensolution with energy E as 

oo 

<r|£)=J>|tr (« 2 )i)(tr (o»)'"|B) = 

where (^(-E 1 ) is some function of the variable E which was not determined by the 
recursion relation eq. (J7j) . Choosing a(n) — 2~ 2j_1 and setting 2E = k 2 we can 
use a known formula for the sum of products of associated Laguerre polynomials 
of the same index [19] , 

Next, by exploiting some properties of the Bessel functions [15], and setting 
<f>(E) — e , we can transform the above result into 

2 2 V 7r kr 

which is, up to a multiplicative factor, the Claudson-Halpern solution of the 
SU{2) model. 

4.2 SU(3) model 
4.2.1 Recursion relation 

We now decompose \E) in the Fock basis of the SU(3) model 

\e) = y, a ^ E ) tr ( af2 y tr (° t3 )' c i )- 

2j+3k<N cut 
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For \E) being an eigenstate, a^k must obey the following recursion relation [TB] . 

Oi-i,* - (2j + 3fc + 4 - 2£)oj, fc + {j + + 3fe + i)a j+ i jk 

+ -(k + l){k + 2)aj- 2 ,k+2=0. (8) 

Notice that the first three terms are diagonal in the k index and are similar 
to the recurrence relation eq.0. The last term in eq.® mixes the coefficients 
with different values of the k index. However, coefficients a,j,k with even and 
odd k remain not related. Therefore, we can solve separately for the amplitudes 
dj,2k and Oj-,2fc+i- 



4.2.2 Finite cut-off solutions 

For reasons of clarity we consider here only the situation when the cut-off N cut 
is even and the solutions contain an even number, 2m, of cubic bricks (a^ 3 )- 
The derivation of solutions with an odd number of cubic bricks is similar. 

It can be shown [18] that the solutions to eq.© can be classified into sev- 
eral separate sets. A solution belongs to the set f m if aj t k = 0, k > 2m and 
a>j,k 7^ 0, k < 2to. In words, this means that the eigenstate can be decomposed 
into basis states containing at most 2m cubic bricks, fo is the simplest set of so- 
lutions, for which only a^o are nonzero, i.e. is only composed of bilinear bricks. 
Each set of solutions has its separate quantization condition for the possible 
values of the E parameter. 

We will now give the general form of solutions belonging to the set f m . 
The set f m contains d m = ^(N C ut ~ 6m)) + 1 solutions with E such that 
L e d ™ +3 (2E) = 0. They can be written s$ [32] 

\ E ) = E L f n+3 ( 2i? )(^ 2TO ) + E r£ ( TO 'P)ij+ 3 ^ 2m - 2 p) 

3=0 P=l 



with 



1 t{2t-\) 



T E (m,p) = Y\--- ^ 

v .w 11 3 2m+ i 



; _ _(2m + l)2-(2f-l)2 

In order to illustrate these formulas, let us present explicitly the solutions from 
the sets fo and f\. They read, respectively, 



\E) =a ^i, 3 (2i?)b-,0), 

3=0 

\E)=a^L](E){\ ]l 2)~^\ ]+ i 1 Q) 



3=0 



3 We adopt a simplified notation in which (tr (a^ 2 )) (tr (at 3 )") |0) = \j,k). 
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The new feature of the />o solutions is the degeneracy, which appears because 
several states can contain the same total number of quanta, i.e. the equation 2j+ 
3k — ub can have several solutions. Particularly, the degeneracy of the states 
containing 6 quanta, namely |2,0) and |0,3) is responsible for the structure of 
the solutions from the set j\ . Figure 2] demonstrates graphically the structure 
of these solutions. 

The complete solution to the eigenvalue problem, i.e. the set of all eigenstates 
{ | E) } is given by the union 

[iN cut \ U( W «««- 3 )J 

{\E}} - (J /*» U U 9m, (9) 

m— m— 

where g m are the corresponding sets of eigenstates with an odd number of cubic 
bricks, and d' m their multiplicities. The spectrum, i.e. the set of all values of 
the E parameter, {-B}, for which a nonzero eigenstate exists, can be written as 

{E}= (J {L^+ 3 (2£)=0} U (J {l%T +6 {2E) =0}. (10) 

m— m=0 

4.2.3 Infinite cut-off solutions 

Eventually, one can show that the solutions retain their structure in the infinite 
cut-off limit. The proof relies on the observation that the mixing coefficients 
T E (m,p) do not depend on N cut . Therefore, in order to obtain the exact solu- 
tion, the sum over j can be safely extended to infinity, 

oo m 

\E) = ]T L 6 J m+3 (2E){\j, 2m) + £ T E {m,p)\ 3 + 3p, 2m - 2p>) , 

i=o p=l 

where E can be now any real number. Notice that in the limit N cut — > oo the 
number of separate sets of solutions f m will become also infinite. 

5 Conclusions 

In this paper we have described the cut Fock space approach to D = 2, su- 
persymmetric Yang-Mills quantum mechanics. We have briefly presented the 
numerical algorithm as well as numerical results, namely the spectra of the 
SYMQM system with the SU(3) gauge group in the bosonic and n F = 2 sec- 
tors. Subsequently, we showed that the cut Fock space approach is also a conve- 
nient framework for analytic calculations. We have derived exact solutions for 
the SYMQM system with the SU (2) gauge group and compared them with the 
original solutions of Halpern and Claudson. Then, we have applied the method 
to the SU(3) SYMQM system and obtained the spectra and eigenstates in the 
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1 2 3 4 5 6 7 

j 

Figure 4: The structure of the solutions of the recursion relation eq. ([8]). Each 
dot represents a coefficient cij.k with appropriate values of the j and k indices. 
The oblique, straight line denotes a cut-off with a fixed number of quanta, here 
Ncut = 15- The coefficients lying below and on this line are included in the 
Fock basis, whereas coefficients lying outside are not. By increasing the cut-off 
we push this line to the right and include more states into the cut Fock basis. 
The remaining lines represent the sets of amplitudes of particular solutions 
of the recursion relation. The lowest, dotted line corresponds to the solution 
involving only quadratic bricks i.e. the solution from the set /o- The dashed 
triple represents a solution from the set f\ . The two horizontal parts of the triple 
denote the amplitudes Oj-,2 and aj.Q. The mixing of these amplitudes starts at 
the number of quanta equal to 6, i.e. both the amplitudes ao,2 and 0,3.0 contain 
6 quanta. 



bosonic sector. Hence, for a given cut-off, one can explain analytically all fea- 
tures of figure [21 

Our analytic results can be extended in several directions. First of all, it 
is possible to obtain recursive relations and solve them in all fcrmionic sectors 
of the model with ST/ (3) gauge group. Obviously, an exact, complete solution 
of this model enables one to calculate the Witten's index [15] . Second, the 
knowledge of exact solutions can be helpful for the investigation of systems 
with interactions. One can use the free solutions at finite cut-off as a starting 
point of perturbative expansion. In both these problems the generalization 
to other SU(N) gauge groups can be achieved. Third, the exact form of the 
solutions in the bosonic as well as fermionic sectors enables one to study their 
large- N limit, which is an important point in the investigations of the SYMQM 



14 



systems. Last but not least, the method can be extended to higher dimensional 
systems, its application to D — 4, SYMQM with SU(2) gauge group is now 
being investigated. 
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